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Abstract 



We investigate a flavored washout effect due to the decay of the lightest right-handed neu- 
trino, assuming that there is non- vanishing initial lcpton asymmetry and the decay of the lightest 
(yr^ I right-handed neutrinos gives negligible contribution to the asymmetry. We figure out general 

00 ■ features of the washout effect. It is shown that there is a novel parameter region where an effect 

that is negligible in most cases plays a critical role and a sizable lepton asymmetry can survive 
against the washout process even in a strong washout region. 
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1 Introduction 



Leptogenesis [T] is a simple mechanism to generate baryon number asymmetry of the Universe. 
The idea is that a lepton asymmetry produced at a high temperature is converted to the baryon 
asymmetry through the sphaleron interactions [2] which conserve B — L number but break B + L 
number. A simplest version of the leptogenesis is based on the seesaw mechanism [3] which can also 
explain the smallness of the neutrino masses by introducing heavy right-handed neutrinos (RHN) 
to the standard model. In the seesaw models, the CP-violating and out-of-equilibrium decay of the 
RHN can produce the lepton asymmetry. 

The seesaw mechanism is also easy to be implemented in supersymmetric and/or grand unified 
theories (GUTs) which are most attractive candidates for the physics beyond the standard model. 
In such a class of model, the successful leptogenesis is considered as a mechanism to generate the 
baryon asymmetry of the Universe. However, in many models especially in GUT models, the mass 
(M\ ) of the lightest RHN (N\ ) is too small to generate the observed baryon asymmetry by the Ni 
decay and the asymmetry should be produced through another mechanism. 

Recently it is pointed out that flavor effects give significant contributions to the leptogenesis 
[H E [7J [8] . One of the interesting phenomena in the flavored leptogenesis is that the primordial 
lepton asymmetry generated by the decay of the second lightest RHN (A^), of the inflaton or so 
can remain against the washout by the lightest one[6j[9]. This is an interesting possibility to give 
enough baryon asymmetry even when the mass of the lightest RHN is too small. In such a scenario, 
the study of the washout effect by the lightest RHN is very important. 

In this paper, we study the detail of this flavored washout effect due to the lightest RHN and 
we point that there is a novel parameter region where an effect that is negligible in most cases 
plays an important role. This effect is due to off-diagonal elements of the so-called A-matrix, and 
thus unique in the flavored leptogenesis. Most recently, the effects of the off-diagonal elements are 
studied numerically in the context that N\ decay generates the lepton asymmetry [10], Here, we 
adapt the analysis on the washout effect to the case where the asymmetry produced by the N2 
decay dominates the baryon asymmetry of the Universe, and show that a sizable lepton asymmetry 
can remain against the washout process in a different way from those studied in Refs.[6j[9]. 

In the section 2, we investigate the washout effect by the N\ decay, assuming that there is 
primordial lepton asymmetry before the decay becomes relevant. In the section 3, we will show 
examples that the primordial asymmetry is generated by the decay of the second lightest RHN. 
The section 4 is devoted to summary and discussion. 

2 Flavor dependence of the Washout Effect 

2.1 The seesaw mechanism and the leptogenesis 

In the seesaw mechanism, RHNs are introduced to the standard model, 

C = C SM + Yfjh'lfNj - ^-NiNf (i = 1, 2, 3 and / = e, //, r) , (1) 

with Ni, lj, and h being RHNs, lepton doublets, and Higgs doublet respectively. Here we take the 
basis where the Yukawa matrix for charged leptons and the mass matrix of RHNs are diagonalized. 
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Integrating out the heavy RHNs and giving the VEV to Higgs, one can obtain the neutrino masses, 
rrii, and Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix [IT], U, as 

ETdiag(mi, m 2 , m 3 )U = u 2 Fdiag(Mf 1 , Mf 1 , M 3 _1 )F T . (2) 

where v = 174GeV is the Higgs VEV. Supposing the reheating temperature after the inflation is 
enough high, the RHNs are produced through the interaction with the doublet leptons and the 
Higgs fields. When the temperature decreases down to the mass of RHN, the production becomes 
inefficient and RHNs decay away. This out-of-equilibrium decay of the RHNs generates B — L 
asymmetry which is proportional to the CP violation in the decay, e, defined as 

j _ ^Ni^ifh - iy^s 

e i — 7 c- W 

l2f {^Ni^lfh + ^Ni^lfh ) 

This asymmetry is converted to the baryon asymmetry through the electroweak sphaleron[2] pro- 
cesses. 



2.2 Boltzmann Equation 

In order to evaluate the baryon asymmetry of the Universe, the Boltzmann equation is used. In this 
analysis, for simplicity, we omit the scattering effects which are considered to be subdominant. The 
decays and inverse decays, N\ <-> Ifh, Ifh, are considered with rate 7^,. With this simplification, 
the evolution of the asymmetry of Af = B/3 — Lj after the decoupling of the second lightest RHN 
is described by the following set of Boltzmann equations [8]: 



Yn, = 
dz 

dyA f z 
dz ~ ~ sH(M\ 




f [ Y m n \ , r D ( vif , vh 



(4) 



(5) 



where z = M\/T and jjj = YlfiD m ^ ne P arame t ers Yx an d Y^ indicate the number density 
of the particle X divided by the entropy density s = 27r 2 g£ ff T 3 /45 and its value in equilibrium 
respectively, and yx = Yx — Y x . The parameter g^ s ~ = 106.75 is the total effective number 
of the degrees of freedom (DOF) at the temperature around M\. With these definitions, we have 

eq 3 45C(3)5f A r 1 j Y£ 1 g N 2 f 1 for fermion 

^-i-S^T'^W. YZ^, ~ 3/4 for boson < 6 > 

where 1S the Riemann's zeta function and K v (x) is the modified Bessel function. Here, gx is 
(not effective) number of DOF of the particle X, for example gi f = gj = 2 and g^ = 2. 

After neglecting the finite temperature effects such as the thermal masses and running of the 
couplings (for these effects, see Ref.[II]) for simplicity, one can obtain 7^ and the Hubble parameter 
H(z) as 



TD = 8Yfi^%T D , H{T) = \r-^ — : (7) 
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where T D = (yty) n Af 1 /(8vr) is the total decay width of Nt and M pl = 1.22 x 10 19 GeV is the 
Planck scale. Now, let 
mass parameter" m* as 



Planck scale. Now, let us define the "washout mass parameter" rhf and "equilibrium neutrino 



mi= ^r> m *= r D = V^o~^r = L07meV ' (8) 

where rhi = Ylf^i- The partial decay width to a flavor / is written as T f D = m{Mf/(8TTv 2 ) 

and the total decay width is given by the sum of them Tp = ^ . Eventually, the Boltzmann 
equations are written as 



dz K 2 {z) m 

s 



dyAf KAz 

-z 



dz K 2 {z) 



m* V JVl / 4 m* \ g lf 2 g h 



(10) 



The coefficient of y^ is different from that in Ref.[8] by a factor 3/4 which comes from the relative 
factor of the number density in the equilibrium of fermion/boson ©. Notice that, however, in the 
derivation of the Boltzmann equations (jU) and ([5]), an approximation / = 1/ (exp ((E — fi)/T) ± 1) ~ 
exp (— (E — n)/T) is made. Within this approximation, the relative factor 3/4 (and the factor 2 
in Eq. (|15p ) disappears. In fact the right hand side of Eq.© is originally written in terms of not 
yx/Y x q but the chemical potentials of the particle X, fix- In literatures, these chemical potentials 
are replaced as Eq.© using the relation between fix and yx/Y^ 1 with the approximation. If one 
would not use this replacement, a factor 1/2 appears instead of 3/4. This difference of the factor 
does not affect the results a lot in many cases and the term of y^ itself is often neglected. In our 
analysis, the contribution can affect the result significantly. In the following we take basically the 
factor 3/4 for illustration. It is straightforward to make analyses with a different factor. 

An important point is that y& f is invariant under the standard evolution of the Universe and 
related to the present baryon asymmetry as ys = 12/37 x V va f [15^1 at the weak scale due to 
the electroweak sphaleron process. Thus, we define "baryon asymmetry" by multiplying the factor 
12/37 on y& f even at a higher temperature. Because this value is proportional to the present baryon 
to photon ratio as ys = g , o ff vr 4 /(45C(3)) x t]b = 7.04ryB where g^ = 43/11 is the present total 
effective number of DOF, successful leptogenesis scenario should predict the "baryon asymmetry" 
y°bs = a87 ± o 03 x 10" 10 which comes from the observable r/° bs = 6.1 ± 0.2 x 10~ 10 [T7]. 

As mentioned in the introduction, in some models, the CP violation t\ is too small to produce 
enough lepton asymmetry. In this case we can neglect the source term. Then, the evolution of y& f 
is controlled by one equation as 

dz 4 m* V gi f 2 g h I 



Note that we assume that the asymmetries of the lepton doublet, y^, and of the Higgs fields, yn, 
are much smaller than 1 and thus neglect the higher terms because it is proportional to the B — L 
asymmetry, as shown below. These relations are forced by the fast (spectator) processes, such as 

1 The factor 12/37 can be somewhat different, for instance 28/79[l6]. depending on the timing of the freeze out 
of the electroweak sphaleron, but in any case, the value is approximately equal to 1/3. 
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the sphaleron process, and depend on the temperature. For example at a high temperature, only 
the interactions mediated by the gauge and the top Yukawa coupling are in the thermal equilibrium, 
while at a lower temperature weaker interactions come in it. For instance, let us concentrate on the 
range of the temperature where the interactions mediated by all the second and third generational 
Yukawa couplings are in the equilibrium but the first generational ones are not. This range is likely 
the one in which we are interested, namely T ~ M\ < 10 9 GeV. In this range, the weak sphaleron 
and the QCD sphaleron [18] are considered to occur fast enough. 

These fast interactions make the following relations hold among the chemical potentials : 

Mft ~~ fJ-Uj + HH = 0, i = 1, 2, 3 (due to the CKM mixing) 

- fJ-dj -fJ>H = 0, j = 2, 3 

~ Me, - m = o, 

Z^i ( 3 M<ft + Wi) = 0, EW sphaleron 

J2i (2/t% - [i Ui - fJLdi) = 0, QCD sphaleron. 

In addition to these relations, we impose the charge neutrality of the Universe and assume the 
vanishing asymmetries for the right handed leptons (and quarks if the QCD sphaleron is not 
considered) of the first generation: 

+ 2 /*m ~ Vck ~ Mi ~ fJ> ej ) + = 0, (12) 

i 

Hei = 0, (13) 

/Aii = Mdi( = ^ no QCD sphaleron) (14) 

Notice that in the Eq. (|12p . the factor 2 in front of /J>h comes from the relative factor in the relation 
between the asymmetry density and the chemical potential for massless particles : 

dxUXrWi f V 2 for fermion . , 

VX = — T { 1 for boson " (15) 

Taking care of this factor 2, we get similar relations among tliG asymmetries by replacing /ifermion — ^ 
2/fcrmion/fi'fcrmion, tyh/ 1 9h and (\ x 3 x (2/J, gi + + fJ, d J) /3 - (2(i lf + /i e/ ) ->■ y Af - By 
solving these relations, we find the expression of yi f /gi f and r/h/dh i n terms of y Ai as 

— = Qff'yA f ,, -.— = -rC hf y Af , (16) 
9l t f 4g h 4 J 1 

with 

/-109/253 25/506 25/506 \ /-53/506\ 

Ci = 29/1012 -493/1518 13/1518 , C h = -37/253 (17) 

\ 29/1012 13/1518 -493/1518/ \-37/253/ 

if we do not consider the QCD sphaleron and with 

/-151/358 10/179 10/179 \ /-37/358\ 
Ci = 25/716 -172/537 7/537 , C h = -26/179 (18) 
\ 25/716 7/537 -172/537/ \-26/179/ 

if we take into account it. In the following, we examine only the latter case because there are no 
qualitative difference. 
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Figure 1: The new variable z' as a function of z. 



From these expressions, we have the following Boltzmann equation: 



dyA f 



Q ~ 7 

Kl ( z )_L A y (19) 

A m . J J J 



dz 4 m* 

with 

/715/716 19/179 19/179 \ / 1.00 0.11 0.1l\ 

Affi = 61/716 461/537 103/537 = 0.085 0.86 0.19 , (20) 

\ 61/716 103/537 461/527/ \0.085 0.19 0.86/ 

where the summation over /' = e, fi, r should be understood. In order to analyse this equation, it 
is convenient to change the variable from z to z' that satisfy ^f- = z 3 Ki(z)/4 so that 

OyAj fh{ 

^ = -t Aff ' y ^- (21) 

The range of z' is from z'(z = 0) = to z'^ = z'(z = oo) = 37r/8 = 1.18. The relation between 
them is shown in the FigJll This figure shows the washout occurs mostly in the temperature range 
Mi/10 < T < 10Mi. 

For comparison, the Boltzmann equation for the usual one-flavor approximation, which is in 
reality valid only when the temperature is high enough that even the processes mediated by the 
tau Yukawa coupling are out of equilibrium, is given aa^| 

dy A fhi 



-VA, (22) 
oz m* 

where Ua = Y1 VAt an d rh\ = ™i- 



2 Here we neglect the Higgs contribution in order to compare our analysis with those in literatures, though it is 
considered in the next section. 
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2.3 Solutions 



Roughly speaking, the matrix Affi is close to a diagonal one. And thus, we can find an approximate 
solution by a perturbation with respect to rather small off-diagonal elements. Neglecting the off- 
diagonal elements, each initial asymmetry is exponentially washed out. The evolution of the 
asymmetry is given by 

J 



y%]( Z ') = ex P (-^A ff z>)yl f . (23) 



When the off-diagonal elements are switched on, the asymmetry follows 



(i) 

f'H 



r)n y > ~ f ~ f 

dz> ~ m* AffVA f If m* ff VA f 



up to the next leading order of the off-diagonal elements. Inserting the leading order solution 
Eq.pS]), we find 



f & — A f'f - —* A ff v v 7 v J ' 

This expression shows that even if the initial asymmetry of a certain flavor is zero, the asymmetry 
is generated from those of the others although it is suppressed by a small off-diagonal element of 
Aij. For most cases, this order already gives good approximation for the final value of the total 
asymmetry. 

(1) rh{ < m* (/ = e,/i,r) 

In this case, one may expect the washout effect is small and thus the flavor effect can not play an 
important role. However, even in this case, the final total asymmetry can be 2 times larger than 
the one-flavor approximation (See (1) in Fig|2|). 



(2) rh[ > m* and m[ < m* (/ / /') 



In this case, the summation rh\ = ^rh^ is dominated by m( and is larger than m*. It is called 
strong washout region. As shown below, however, if we take account of the flavor effect the 
washout effect is drastically changed. The effect depends strongly on the flavor structure of the 
initial asymmetry. 

In the following, let us take mf < m* 5j rhf (a = fx, r) as a representative example for clarity. 
It is straightforward to apply this analysis to the other cases. We consider two typical sets of the 
initial asymmetries: 

(2a) y° Ae > y% a 

In this case, the washout effect of j/a e is controlled by fh\ which is much smaller than rh\, 
while yA a are generated due to the small off-diagonal elements of A ea , with the opposite sign. 
Because of the large mf, these generated y& a are washed out strongly and can not become 
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comparable with y& e . Thus, in terms of the perturbation, the leading order approximation 
is sufficient. 

Note that even in this case, the washout factor yA{z' 00 )/y A ~ exp (rhf A ee /m*) is quite differ- 
ent (much larger) than the one- flavor approximation, exp (mi /m*) (See (2a) in Figj2]). This 
is the case even when m\ is not so much smaller than the others, due to the exponential 
washout factor [6]. 

(2b) y° Ae < y° A ^ and/or VAt 

The asymmetry y& a decreases rapidly, while the asymmetry y& e produced due to the off- 
diagonal elements is washed out much more slowly. This means that at some point yA e 
becomes dominant. Once it becomes dominant, the following evolution is similar to the one 
in the case (2a). Thus, the washout factor is controlled basically by the small rh\ rather than 
mi or rh\ (See (2b) in Figj2|). Interestingly in this case, the sign of the total B — L asymmetry 
changes through the washout. 

Note that the approximation at the NLO is quite bad for y& a because the secondary con- 
version from j/A e j which is generated by the NNLO effect, is important. Nevertheless, the 
approximation for the total asymmetry is rather good because these are small as in the case 
(2a). 

(3) fh[ > m* (/ = e,n,r) 

In this case, all the asymmetries in each flavors are strongly washed out. Thus, it is hard that the 
observed value remains after the washout, as far as M\ < 10 9 GevH. 

From the above considerations, it is clear that the washout factor is basically controlled by the 
smallest washout mass parameter. This is in great contrast to the non-flavored case, where the 
factor is controlled basically by the largest washout mass parameter. Interestingly, this is also true 
even for the case that the initial asymmetry of the flavor with smallest washout mass parameter is 
tiny. For this case, the effect of the off-diagonal elements, which is usually negligible, is crucial. 

2.4 Fixed point 

From Eq. (f2Tj) . one can obtain a coupled equations for yA f /lJA T (/ = e,//) as 



_d_ fVA f \ m[ 



^A ff ,-A rf ( y -± 
ml '' \yA T 



y A T 



(26) 



This couple of equations have fixed points in the space of (yA e / UA T i UA^ / VA T ) i which is determined 
by solving the equations 



E 

f'=e,n,T 



m\ V y At 



^)=0. (27) 

VA T 



3 If we consider models with M\ ^> 10 9 GeV where the muon Yukawa interaction is out of equilibrium, the 
asymmetry along with the direction in the flavor space that is orthogonal both to the direction determined by Ni 
Yukawa coupling and r direction is free from the washout, even in this case [9]. 
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(2b) {mlm^ml} = {2,15,20}meV , {y° Aa , y° A ^ y° Ar } = {0,0.5,0.5} 

Figure 2: The evolutions of the B/3 — Lf asymmetries. The horizontal line is z'/ z oo' anc ^ the 
vertical line is log 10 |t/a| • In the left figures, black (solid), light blue (broken) and purple (dot- 
dashed) lines show the total B — L asymmetry calculated by the full Boltzmann equation (|21|) . by 
the approximation formula (I25|) and the one-flavor approximation (j22[) . respectively. In the right 
figures, red, green and blue lines respectively show the A e , and A T , and the solid and broken 
lines corresponds (i2~Tj) and (i25|) . respectively. 
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Once the flow of the solution reaches close to a fixed point at z' = z[ , the set of ratios (yA e /va t , 2/A M / VA T ] 
becomes invariant and the Boltzmann equations can be rewritten as 

< 28 > 

This means that the asymmetries of the all flavor are washed out with the universal washout mass 
parameter which corresponds to one of the eigenvalues of matrix Bffi = rh^Atft. There 
are three possible points in the case of three effective flavor numbers. However two of them are 
unstable fixed points and only one point is attractive. The attractive one most likely corresponds 
to the smallest eigenvalue which is smaller than the smallest m{. Around the attractive fixed point, 
the asymmetry in the flavor with the smallest washout mass parameter, VAj dominates the total 
asymmetry. This can be understood as follows. The asymmetry in the flavor with larger washout 
mass parameter, VA f decrease more quickly as discussed in the case (2b) in the Sec. 12.31 When 
UAf becomes the order of A^^UAf , UA f2 evolves similar to VA fl because the transportation from 
UAf becomes to control the evolution. Then the washout of VA f becomes weaker due to the 
transportation from UA f2 ■ 

When the initial condition is too far from the fixed point and fh\ <C m*, the washout term 
decouples from the system before the solution flows into the fixed poinlQ. 

One can see similar phenomena also in the case where only N\ decay produces the asymmetries 
and they are washed out by the N\ (inverse) decay. In fact effects of off-diagonal elements of A 
matrix in such a case are discussed in Ref. [10| and the asymmetry shown there behaves as discussed 
in the above. 



3 Asymmetries by the Second Lightest RHN Decay 

In this section, we investigate the possibility that the initial asymmetries are generated via the 
second lightest RHN decay. 

For simplicity, we restrict ourselves to the case that the mass of the second lightest RHN is 
larger than 10 12 GeV. This case is qualitatively discussed in Ref. [9]. For this mass range, the fast 
interactions that are in the equilibrium when the RHN decays are only the interactions mediated 
by the top Yukawa coupling and QCD sphalerons. This means that the fast interactions can 
not distinguish all the generations of the lepton doublets, and thus two liner combinations of 
the three doublets that do not interact with the RHN are never produced. Namely, only In oc 
Yti^t + Yuzln + Y e 2le are produced. Then, the relevant Boltzmann equations are for one flavor 
system, which is given by Eqs.(|9|) and (I10p by replacing all the index 1 to 2 (including z — > M2/T) 
and suppressing the flavor indexes. With the definitions given in the section 12.21 we have 

The A factor for this case is calculated in a similar way to the discussion in the section 12.21 as 

A = 2Q + \c h = g (31) 



4 Notice that the z' takes the value in the range of < z' < 3n/8 with respect to < 2 < 00. 
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It is possible, of course, that we solve these set of equations numerically to evaluate the B — L 
asymmetry produced by the decay of the second lightest RHN. In this article, however, we use the 
following approximation formula proposed in Ref.[19j, which includes the effects of the scatterings, 
to evaluate the "baryon asymmetry": 

12 

VB ^^e2v(Arh 2 ) (32) 



with 

x 



-1.16\ _1 



In any case, these equations are controlled by the parameters mi and £2 which are determined 
by the mass spectrum of the RHNs Mi and the neutrino Yukawa coupling Yfi. To be more concrete, 
they are given by sums of ([8]) and 

/ll, /7„ t „\ „fM?\ / M} 



'2 



Here both the diagrams of the vertex correction and of the self-energy correction are implemented 
in each function as 

= + V£ (l - (1 + *) In > ( 35 ) 

In this way, fixing Mj and Yfi, all the parameters in the Boltzmann equations are determined, 
and we can calculate the B — L asymmetry i/a just after the second lightest RHN decouples. 
This asymmetry is washed out when the lightest RHN starts decaying. In this period, the r and 
\x Yukawa couplings enter into the equilibrium, and thus the fast interactions distinguish all the 
three flavors. Therefore, we should divide da into ua t , VA^ and yA e , which follow the relation 
VA T '■ 2/A M : VA e = "^2 • 1 ^ l 2 '■ ™2 according to the probabilistic interpretation. This set of 
asymmetries UA f , (f = r) gives the initial condition of the analysis given in the section [2T3l 

The neutrino Yukawa couplings Y should be related to the low energy neutrino parameters 
through the seesaw relation, Eq.©. In order to represent the solution of this relation, we adopt 
the following famous parameterization [20 j . 



Y f i = (U*) fj ^mjRji^/W/v. (37) 

Here U is written as the product of a CKM-like mixing matrix V which includes three mixing angles 
and one CP phas^land a phase matrix with two Majorana phases P = diag (1, exp (2021/2) , exp [ia^\/2)) 
: U = VP and R is a complex orthogonal matrix which can be decomposed as R = e tuJ23X? e luJlsXs e^ 12 ^ 2 
where Aj are Gell-Mann matrices and u>ij are complex parameters. For simplicity, in this article, 
we use the following set of the parameters for the light neutrino sector as rrii = {0, 9, 50} meV, 
{sl2 2 , s23 2 , sl3, 5, (X21 , 031 } = {0.3,0.5,0,0,0,0} for the PMNS matrix, and Majorana masses 
Mi = {10 7 , 10 13 , 10 14 } GeV for the RHNs. 



As a parametrization of V, we adopt the Chau-Keung parametrization (PDG parametrization) [21] 
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(I) (la) 



rh (meV) 



e(l(T 6 ) 





y^(10- 10 ) {0.97,1.26,1.13} {2.53,2.09,4.28} 

ysilO" 10 ) 1.26 1.02 




Table 1: Results of rh, e, Yg and Yb for the three examples in ([38 



As representative examples, let us consider the following sets: 



(I) 
(Ea) 

(2b) 



{^12, ^23, ^13} 
{U>12, UJ23, ^13} 
{^12,^23)^13} 



{30°,i5°,-l°} 

{-88°,(60 + i3)°,3°} 

{(-85 + »4)°, (50 + i20)°, -5.5°} 



(38) 



For instance, for the example (I), we find 



m 




2.04 
2.66 
2.39 




meV, 



-0.27 -0.02 
-95.0 11.57 
84.3 -11.28, 



x 10" 



(39) 



These show {mf,m^,?n[} = {0.68, 0.71, 0.98} meV < m* and {ef , e^*, e^} = {lO" 13 , 10" 11 , 10~ n } 
are negligibly small. Thus, this is an example of the case (1) in the section [2~3l Using the approx- 
imation (|32p . we see the "baryon asymmetry" generated by the N2 decay is Yg = 3.36 x 10 -10 . 
When the temperature decrease to around Mi, Ni starts decaying, and the asymmetry is washed 
out in the way investigated in the last section. As mentioned above, in this period, the fast in- 
teractions distinguish all the flavor, and the asymmetry should be divided as jyg e , y^, y# T j = 



{0.97, 1.26, 1.13} x 10 . After the washout, a total asymmetry ys = 1.26 x 10 remains. 

In a similar way, (Ha) and (lib) are examples of the cases (2a) and (2b), respectively. Their 
results are listed in the Table 



4 Summary and Discussions 

In this article, we investigate the washout effect due to the N± (inverse) decay, assuming non- 
vanishing initial lepton asymmetry and negligible lepton asymmetry production in N\ decay. We 
show that there is a novel parameter region in addition to those studied in Refs.jU [9]. There, 
off-diagonal elements of the ^-matrix, which are often omitted, play a critical role. This region is 
where some of m[ is comparable to or smaller than m*, the others are larger than it, and the initial 
asymmetries on the flavors with small m{ are tiny. In this case, if we would omit the off diagonal 
elements as usual, any initial asymmetries on the flavors with large m{ were strongly washed out. 
In fact, the off diagonal elements transform the asymmetries from those with large m{ to those 
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with small ones. Once transformed, such asymmetries are weakly washed out, and thus a sizable 
total asymmetry may survive. 

For completeness, we examine the possibility that the initial asymmetry is generated by the N2 
decay within the thermal leptogenesis scenario. We show some concrete examples for each class 
discussed in the above analysis. 

Finally, let us make a comment on an ambiguity of the Boltzmann equations, especially on 
the factor in front of y^ in Eq. (|10p . As briefly discussed below the equation, an approximation is 
used in the derivation of the Boltzmann equations ([4]), (|5|), and it brings the ambiguity. Because 
the contribution of this term (y^) is relatively small, as seen from (|17p and (|18|) . this ambiguity 
does not affect results so much, 0(10%). In the case (2b), however, the off diagonal element of 
A-matrix, which is of the same order as the y^ contribution, is critical. In addition, a cancellation 
occurs between y^ contribution and that of yi in the off diagonal element. Thus, the result is largely 
changed due to the ambiguity. In fact, if we take a factor 1/2 instead of 3/4 as a possible choice, 
the final baryon asymmetry is reduced to ys = 0.83 x 10 -10 with the same parameters as (lb) in 
(|38|) . Thus, it is important to make a closer look on the Boltzmann equations before discussing 
this novel effect quantitatively. 
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